Keywords

We present the results of PSRπ, a large astrometric project targeting radio pulsars using the Very Long Baseline Array (VLBA). From our astrometric database of 60 pulsars, we have obtained parallax-based distance measurements for all but 3, with a parallax precision that is typically ∼45 μas and approaches 10 μas in the best cases. Our full sample doubles the number of radio pulsars with a reliable (≳5σ) model-independent distance constraint. Importantly, many of the newly measured pulsars are well outside the solar neighborhood, and so PSRπ brings a near-tenfold increase in the number of pulsars with a reliable model-independent distance at d > 2 kpc. Our results show that both widely used Galactic electron density distribution models contain significant shortcomings, particularly at high Galactic latitudes. When comparing our results to pulsar timing, two of the four millisecond pulsars in our sample exhibit significant discrepancies in their proper motion estimates. With additional VLBI observations that extend our sample and improve the absolute positional accuracy of our reference sources, we will be able to additionally compare pulsar absolute reference positions between VLBI and timing, which will provide a much more sensitive test of the correctness of the solar system ephemerides used for pulsar timing. Finally, we use our large sample to estimate the typical accuracy attainable for differential VLBA astrometry of pulsars, showing that for sufficiently bright targets observed eight times over 18 months, a parallax uncertainty of 4 μas per arcminute of separation between the pulsar and calibrator can be expected.

Magnetars are a class of highly magnetized, slowly rotating neutron stars, only a small fraction of which exhibit radio emission. We propose that the coherent radio curvature emission is generated by net charge fluctuations from a twist-current-carrying bundle (the j-bundle) in the scenario of magnetar-quake. Two-photon pair production is triggered, which requires a threshold voltage not too much higher than 10⁹ V in the current-carrying bundle, and which can be regarded as the "open field lines" of a magnetar. Continued untwisting of the magnetosphere maintains change fluctuations, and hence coherent radio emission, in the progressively shrinking j-bundle, which lasts for years until the radio beam is too small to be detected. The modeled peak flux of radio emission and the flat spectrum are generally consistent with the observations. We show that this time-dependent, conal-beam, radiative model can interpret the variable radio pulsation behaviors and the evolution of the X-ray hot spot of the radio-transient magnetar XTE J1810−197 and the high-B pulsar/anomalous X-ray pulsar PSR J1622−4950. Radio emission with luminosity of $\lesssim {10}^{31}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$ and high-frequency oscillations are expected to be detected for a magnetar after an X-ray outburst. Differences of radio emission between magnetars and ordinary pulsars are discussed.

The second repeating fast radio burst source, FRB 180814.J0422+73, was detected recently by the CHIME collaboration. We use the ten-year Fermi Large Area Telescope archival data to place a flux upper limit in the energy range of 100 MeV−10 GeV at the position of the source, which is ∼1.1 × 10⁻¹¹ erg cm⁻² s⁻¹ for a six-month time bin on average, and ∼2.4 × 10⁻¹² erg cm⁻² s⁻¹ for the entire ten-year time span. For the maximum redshift of z = 0.11, the ten-year upper limit of luminosity is ∼7.3 × 10⁴³ erg s⁻¹. We utilize these upper limits to constrain the fast radio burst (FRB) progenitor and central engine. For the rotation-powered young magnetar model, the upper limits can pose constraints on the allowed parameter space for the initial rotational period and surface magnetic field of the magnetar. We also place significant constraints on the kinetic energy of a relativistic external shock wave, ruling out the possibility that there existed a gamma-ray burst (GRB) beaming toward Earth during the past ten years as the progenitor of the repeater. The case of an off-beam GRB is also constrained if the viewing angle is not much greater than the jet opening angle. All of these constraints are more stringent if FRB 180814.J0422+73 is at a closer distance.

The high sensitivity and wide frequency coverage of the Murchison Widefield Array allow for the measurement of the spectral scaling of the pulsar scattering timescale, α, from a single observation. Here we present three case studies targeted at bright, strongly scattered pulsars J0534+2200 (the Crab pulsar), J0835−4510 (the Vela pulsar), and J0742−2822. We measure the scattering spectral indices to be −3.8 ± 0.2, −4.0 ± 1.5, and −2.5 ± 0.6 for the Crab, Vela, and J0742−2822, respectively. We find that the scattered profiles of both Vela and J0742−2822 are best described by a thin screen model where the Gum Nebula likely contributes most of the observed scattering delay. For the Crab pulsar we see characteristically different pulse shapes compared to higher frequencies, for which none of the scattering screen models we explore are found to be optimal. The presence of a finite inner scale to the turbulence can possibly explain some of the discrepancies.

We study the gamma-ray emission patterns and light curves in dissipative pulsar magnetospheres. We produce the gamma-ray light curves by using the geometric method and the particle trajectory method. For the geometric method, assuming the gamma-ray emission originates in a finite-width layer along the last closed lines, we generate the gamma-ray light curves based on the uniform emissivity along the magnetic field lines in the comoving frame. For the particle trajectory method, we consider the spatial distribution of conductivity σ by assuming a very high conductivity within the light cylinder (LC) and a finite conductivity outside the LC. Assuming that all the γ-ray emission originates in the outer magnetosphere outside the LC, we generate the gamma-ray light curves by computing realistic particle trajectories and Lorentz factors, taking into account both the accelerating electric field and curvature radiation loss. Further, we compare the modeling light curves to the observed light curves at >0.1 GeV energies for the Vela pulsar. Our results show that the magnetosphere with the low σ value is preferred for the geometric method. However, the magnetosphere with a near force-free regime within the LC and a high σ value outside the LC is favored for the particle trajectory method. It is noted that the particle trajectory method uses the parallel electric fields that are self-consistent with the magnetic fields of the magnetosphere. We suggest that the results from the particle trajectory method are better supported on the physical ground.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars (MSPs) without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the MSPs PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821−24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the σ_z measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities σ_z ≈ 3 × 10⁻¹⁴ for PSR B1937+21 and on the order of 10⁻¹² for PSRs B1821−24 and J0218+4232. Within the span of our NICER data we do not yet see the characteristic break point in the slope of σ_z; detection of such a break would indicate that further improvement in the cumulative root-mean-square timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821−24 and PSR B1937+21 on timescales of >2 yr.

We present the earliest X-ray observations of the 2018 outburst of XTE J1810−197, the first outburst since its 2003 discovery as the prototypical transient and radio-emitting anomalous X-ray pulsar (AXP). The Monitor of All-sky X-ray Image (MAXI) detected XTE J1810−197 immediately after a November 20–26 visibility gap, contemporaneous with its reactivation as a radio pulsar, first observed on December 8. On December 13 the Nuclear Spectroscopic Telescope Array (NuSTAR) detected X-ray emission up to at least 30 keV, with a spectrum well-characterized by a blackbody plus power-law model with temperature kT = 0.74 ± 0.02 keV and photon index Γ = 4.4 ± 0.2 or by a two-blackbody model with kT = 0.59 ± 0.04 keV and kT = 1.0 ± 0.1 keV, both including an additional power-law component to account for emission above 10 keV, with Γ_h = −0.2 ± 1.5 and Γ_h = 1.5 ± 0.5, respectively. The latter index is consistent with hard X-ray flux reported for the nontransient magnetars. In the 2–10 keV bandpass, the absorbed flux is 2 × 10⁻¹⁰ erg s⁻¹ cm⁻², a factor of 2 greater than the maximum flux extrapolated for the 2003 outburst. The peak of the sinusoidal X-ray pulse lags the radio pulse by ≈0.13 cycles, consistent with their phase relationship during the 2003 outburst. This suggests a stable geometry in which radio emission originates on magnetic field lines containing currents that heat a spot on the neutron star surface. However, a measured energy-dependent phase shift of the pulsed X-rays suggests that all X-ray emitting regions are not precisely coaligned.

Testing gravity using binary pulsars has become a key contemporary focus. Screened modified gravity is a kind of scalar-tensor theory with a screening mechanism in order to satisfy the tight solar system tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the 10 binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS)–white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NSs and WDs are PSR B1913+16 and PSR J1738+0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

The anomalous X-ray pulsar XTE J1810−197 was the first magnetar found to emit pulsed radio emission. After spending almost a decade in a quiescent, radio-silent state, the magnetar was reported to have undergone a radio outburst in 2018 December. We observed radio pulsations from XTE J1810−197 during this early phase of its radio revival using the Ultra-Wideband Low receiver system of the Parkes radio telescope, obtaining wideband (704–4032 MHz) polarization pulse profiles, single pulses, and flux density measurements. Dramatic changes in polarization and rapid variations of the position angle of linear polarization across the main pulse and in time have been observed. The pulse profile exhibits similar structures throughout our three observations (over a week timescale), displaying a small amount of profile evolution in terms of polarization and pulse width across the wideband. We measured a flat radio spectrum across the band with a positive spectral index, in addition to small levels of flux and spectral index variability across our observing span. The observed wideband polarization properties are significantly different compared to those taken after the 2003 outburst, and therefore provide new information about the origin of radio emission.

An accreting millisecond X-ray pulsar, IGR J17591−2342 was discovered in 2018 August in scans of the Galactic bulge and center by the International Gamma-Ray Astrophysics Laboratory X-ray and gamma-ray observatory. It exhibited an unusual outburst profile with multiple peaks in the X-ray, as observed by several X-ray satellites over 3 months. Here we present observations of this source performed in the X-ray/gamma-ray and near-infrared domains and focus on a simultaneous observation performed with the Chandra High Energy Transmission Gratings Spectrometer (HETGS) and the Neutron Star Interior Composition Explorer (NICER). The HETGS provides high-resolution spectra of the Si edge region that yield clues as to the source's distance and reveal evidence (at 99.999% significance) of an outflow with a velocity of 2800 km s⁻¹. We demonstrate good agreement between the NICER and HETGS continua, provided that one properly accounts for the differing manners in which these instruments view the dust-scattering halo in the source's foreground. Unusually, we find a possible set of Ca lines in the HETGS spectra (with significances ranging from 97.0% to 99.7%). We hypothesize that IGR J17591−2342 is a neutron star low-mass X-ray binary at the distance of the Galactic bulge or beyond that may have formed from the collapse of a white dwarf system in a rare, calcium-rich Type Ib supernova explosion.